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Raman-assisted hopping may be used to realize the anyon Hubbard model in one-dimensional 
optical lattices. We propose a feasible scenario that signihcantly improves the proposal of [T. 
Keilmann et ah, Nature Commun. 2, 361 (2011)], allowing as well for an exact realization of 
the two-body hard-core constraint, and for controllable effective interactions without the need of 
Feshbach resonances. We show that the combination of anyonic statistics and two-body hard-core 
constraint leads to a rich ground state physics, including Mott insulators with attractive interactions, 
pair superfluids, dimer phases, and multicritical points. Moreover, the anyonic statistics results 
in a novel two-component superffuid of holon and doublon dimers, characterized by a large but 
hnite compressibility and a multipeaked momentum distribution, which may be easily revealed 
experimentally. 
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Particles are classified as bosons or fermions depending 
on whether their wave function is symmetric or antisym¬ 
metric under exchange. Other types of quantum statis¬ 
tics are, however, possible in lower dimensions. Remark¬ 
ably, 2D systems allow for the existence of anyons, i.e. 
particles with fractional statistics interpolating between 
bosons and fermions [1-3]. Anyons play a fundamental 
role in key areas of modern physics, as fractional quan¬ 
tum Hall effect [4-7] and topological quantum comput¬ 
ing [8]. Fractional statistics is, however, not exclusive of 
2D systems [9]. In particular, ID anyons have attracted a 
large deal of interest [10-23], although the experimental 
realization of a ID anyon gas is still lacking. 

Ultracold atoms offer extraordinary possibilities for the 
analysis of interesting many-body systems [24]. In par¬ 
ticular, several ideas have been proposed for the creation 
and manipulation of anyons in cold gases [25-29]. Par¬ 
ticularly interesting is the recent proposal for the realiza¬ 
tion of the anyon Hubbard model (AHM) using Raman- 
assisted hopping in ID optical lattices [30]. In this pro¬ 
posal the anyonic statistics may be controlled at will, 
opening the possibility for the observation of statisti¬ 
cally induced quantum phase transitions [30], asymmet¬ 
ric momentum distributions [21], and intriguing particle 
dynamics in the lattice [20, 22, 23]. 

In this Letter we first discuss a scheme for realizing 
the AHM that, although following Ref. [30], solves cru¬ 
cial drawbacks that would render the original proposal, 
in general, unfeasible. The scheme also allows for con¬ 
trollable effective interactions without the need of Fes¬ 
hbach resonances, and for an exact two-body hard-core 
constraint (2BHCC), contrary to the approximate real¬ 
ization of 2BHCC resulting from Zeno projection due 
to large three-body loss rates [31], where non-negligible 
losses are typically present. Neither the controllable in¬ 
teractions nor the inherent 2BHCC were considered in 
Ref. [30], which focused exclusively on statistically in¬ 
duced superffuid (SF) to Mott insulator (MI) transi¬ 
tions. We show that the interplay of anyonic statis- 
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FIG. 1: (a) Raman scheme proposed for the realization of the 
AHM; (b) Raman assisted hops (i)-(iv) discussed in the text. 

tics, 2BHGG, and controllable interactions, results in a 
far richer physics for the AHM, including pair super¬ 
fluid (PSF), a dimer (D) phase, and an exotic partially 
paired (PP) phase. The latter constitutes, to the best of 
our knowledge, a novel two-component superffuid charac¬ 
terized by a large, but finite, compressibility and a mul¬ 
tipeaked momentum distribution, which may be readily 
revealed in time-of-flight (TOF) experiments. 

Anyon-Huhhard model- We briefly describe our pro¬ 
posal for the realization of the AHM [32]. We consider 
atoms, bosons or fermions, with states \A) and |5), in 
a deep tilted spin-independent ID lattice with no direct 
hopping. For the specific case of ^^Rb, the choice would 
be \A) = \F = l^mp = —1) and \B) = \F = 2, mi? = 
—2). Both states are coupled far from resonance by Di 
lasers (with linear polarization) and 1/2,3 (with cir¬ 
cular cr_ polarization) (Fig. 1(a)). Due to selection rules, 
\B) is just affected by lasers Li^ 4 . In contrast, both 1/2,3 
and 1/1^4 couple with |A), but the coupling with Li ^4 can 
be made much smaller than that of 1/2,3, and does not 
affect the main conclusions of this Letter [33]. We hence 
assume below that \A) is just affected by Li^ 2 - These 
lasers, with frequencies c^j=i,...,4, induce four Raman- 
assisted hops (Fig. 1(b)): (i) (A,0) ^ (0,A), given by 
7^2,3, such that UJ 2 — 0 J 3 = —A (with A the lattice tilting); 
(ii) {A, A) { 0 ,AB), assisted by 1/2,4 with uj 2 — oj 4 = 
-A + Uab + U; (iii) {AB, 0 ) {A, A), given by 1/^3 
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FIG. 2: Phase diagram of the AHM: (a) dilute limit (p ^ 0), (b) incommensurate intermediate fillings (here p = 1/2) and (c) 
unit Filing (p = 1). The PSF-SF transition lines (boxes or solid line in (a)), extracted from the crossing of two and one particle 
excitations coincides with the onset of parity order. At incommensurate fillings for small 0 the system exhibits a first order 
transition from SF to PSF phase through a forbidden region (FR, dotted area). For (^/tt > 0.6 we find the PP phase. For 
unit filling the Ising-type PSF-SF transition and the Kosterlitz-Thouless SF-MI transition lines cross around ^/tt 0.8 in a 
multicritical point. For larger values of (j) the system exhibits a dimerized gapped phase D. (d) Grand-canonical phase diagram 
and lines of constant density p for = 0 as function of p and 0. Solid lines denote the MI and PP phase, the rest is SF. For 
p/t < —2 {p/t > 4) the vacuum p = 0 (fully occupied phase, p = 2) is realized. 


with cji — Co’s = — A — Uab — U ; (iv) {AB^ A) (A, AB ), 
assisted by such that ui —004 = —A. The notation 
(Vl^Vr) denotes the state of neighboring sites L and R: 
Vl,r = 0 (vacuum), A (single A occupation) or AB (an A 
and a B particle); Uab characterizes the AB interaction. 
The frequency differences compensate the tilting, avoid¬ 
ing Bloch oscillations; moreover, the detuning U acts as a 
controllable effective on-site interaction, as shown below. 

We assume the width of the Raman resonances, VF, 
is small-enough such that each Raman process may be 
addressed independently. If the atoms are fermions, the 
only undesired process (A, 0) ^ (0, B) (where B denotes 
a site with a single B atom) may be avoided if Uab ^ 
W. For bosons, {A, A) (0,AA) and {AB,A) 

(B, A A) (with A A a site with two A atoms) must be also 
avoided. This demands Uaa — Uab^ Uaa^ Uab ^ 
where Uaa characterizes AA interactions [32]. 

The laser Lj is characterized by a Rabi frequency 
with (j)i = —0, and (pj^i = 0. We impose 

|t^l||t^4| _ |02||t^3| _ _ |02||ri4| _ q 2 87p-U. 

4 ~ 3 “ 2v^ “ 2V^ “ 

for other species Clebsch-Gordan coefficients may dif¬ 
fer [32]). \Qj\ is chosen such that they mimic bosonic 
enhancement, absent otherwise due to the distinguisha- 
bility between A and B; e = 1 (—1) for bosons (fermions). 
Under these conditions the system reduces to an effec¬ 
tive single-component ID Bose-Hubbard model [32], with 
on-site Fock states {|0), |1) = |A), |2) = |A5)}, and an 
occupation-dependent Peierls phase [30, 32]: 

^ + H-C-) + I - 1) (1) 

j j 

where U is an effective coupling constant whose sign and 
magnitude may be tailored without the need of Fesh- 
bach resonances by properly choosing the laser frequen¬ 
cies; t /5)/2/\ is the effective hopping rate [34], 


with J the lattice hopping in the absence of tilting, and 
S ^ Uab, a, fl the detuning from the single-photon tran¬ 
sitions; bj (bj) are creation (annihilation) operators at 

site j, and rij = b^jby Introducing aj = 

one obtains the anyon-Hubbard model (AHM) [30]: 

H = -iEA«i+i+H-c.) + ^ E^iK- -1) (2) 

j j 


where satisfy anyonic commutations [45]: OLjOL\ — 

g-i(/)sgn(j-k)^t^^ _ ajak — = 0. 

None of the Raman-assisted hoppings leads to triple 
occupation, resulting in an effective 2BHCC, (5^)^ = 0. 
This inherent constraint turns out to be crucial. It pre¬ 
vents collapse for a large U/t < 0, leading to a rich 
physics of quantum phases (see Fig. 2) [46]. 

Note that although hops (i) and (iv) share the same 
energy difference —A, they may be independently ad¬ 
dressed, even if J ^ Uab,^- This point is crucial, 
solving the major drawback of Ref. [30] in which only 
one ground state component was considered, and hence 
hops (i) and (iv) cannot be independently addressed un¬ 
less 6 Uab, which for typical experimental param¬ 
eters would lead to large heating. An exception may 
be given by the use of clock transitions in alkaline-earth 
atoms [47]); however, the adiabatic elimination of the ex¬ 
cited state demands U ^ (5, and hence t ^ 6 Uab, A, 
would lead to very long experimental time scales [48]. 

Effective repulsion.- The dilute limit (lattice filling 
p <C 1) is best understood from the analysis of the 
two-particle problem [49], characterized by the scatter¬ 
ing length (in units of the lattice spacing) 


— (1 + cos (f) 

4(1 — cos (j)) + 2U/t 


(3) 


By comparison to a ID gas of bosons with mass m one ex¬ 
tracts an effective interaction strength g = —2j{am) [49]. 
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FIG. 3: Phase transitions for 0 = tt. (a) Equation of state, density p as function of chemical potential p (calculated for L = 48 
sites, U — The inset shows a detailed view on the PP region for L — 240 sites, (b) Estimation of the central charge for 
L = 24 (dashed line) and L = 48 (solid line) sites and periodic boundary conditions (p = 1/2). The inset shows fits to the 
block-entanglement entropy SvN for (from bottom to top) U = 2,—2,—1.0. (c) Sum of local single- (Na) and two-particle 
correlations (Nd) for L = 60 sites and p = 1/2, see text. 


Thus the AHM behaves as a repulsively interacting sys¬ 
tem for small p and U ^0, since even in that limit g 
remains finite and positive for any 0 > 0, approaching 
Tonks limit as 0 ^ tt. This is consistent with a weak 
coupling analysis [50], in which the Luttinger liquid (LL) 
parameter K = 7r/(0^ + ^)^/^, showing that <p has qual¬ 
itatively the same effect as a repulsive [/ > 0. In Fig. 2(a) 
we depict the phases for p ^ 0 [51]: a SF phase, and a 
PSF of bound pairs at <0. Due to the effective repul¬ 
sion the SF extends into the < 0 region for 0 > 0. 

Incommensurate fillings.- Bose models (0 = 0) with 
2BHCC have been recently studied in different scenar¬ 
ios [31, 52, 53]. A SF to MI transition results for p = 1 
and large-enough U > 0; for < 0 a transition to a PSF 
occurs, which at incommensurate p may become first or¬ 
der [54]. This forbidden region (FR) occupies the shaded 
area in Fig. 2(b). At fixed p the ground state is charac¬ 
terized by the formation of a macroscopically bound state 
with a phase separation of doublons and holes. The FR 
shrinks, however, fast with increasing 0. 

For (/) > O.Gtt a novel gapless phase, henceforth called 
partially paired (PP) phase, appears between the PSF 
and SF phases. As depicted in Fig. 3(a), a sharp kink 
in the p{p) curve (at p —t) accompanies the onset of 
the PP phase. Whereas SF (PSF) is characterized by 
steps of 1 (2) particles in the p{p) curve, the PP phase 
exhibits a complex sequence of steps of 2 and 1 particles. 
Even more relevant is the behavior of the central charge c, 
which we calculated from the conformal expression of the 
von Neumann entropy, Svn,l{1) = f sin (J/)] + 7, 
for a subsystem of length I in a system of L sites, with 
7 a constant [55, 56]. As shown in Fig. 3(b), contrary to 
the SF and PSF phases (with c = 1), the PP phase is a 
two-component gapless phase characterized by c = 2. 

Nature of the PP phase.- The features of the PP 
phase may be understood from the following simpli¬ 
fied picture. Due to the 2BHCC the hopping term in 
Eq. (1) is Ji6]&j+i + J 2 b]njbj+i, with Ji/t = 1 and 


J 2 /t = — 1 + For (p ^ TT correlated hopping dom¬ 
inates, J 2 = —2Ji, and doublons, or more precisely 
dimers, \d)j = a \11)\20) gain a large bind¬ 
ing energy — \/2J2. However, dimer hopping, Jil^/2, 
is reduced compared to single particles. For a certain 
range of small U and p it is energetically favorable to 
occupy both doublon dimers \d) ■ and atomic dimers, 

\a). = P jp+v Neglecting interactions 

between these two quasiparticles, reasonable for small 
doublon and atom densities, pd and pa, with paPpd ^ 1, 
one arrives at an effective model with two independent 
hardcore bosons. Had = Ha P Hd, with (for U = 0) 

Ha = Ji + H.C.) - Ji (4) 

j 3 

H, = (5) 

^ 3 3 

where aj (at) and dj (dj) are annihilation (creation) op¬ 
erators of hardcore atom and doublon dimers on sites 
+ 1, and pa P 2pd = p. In spite of being a gross 
oversimplification. Had captures the qualitative features 
of the PP phase. For U ^ 0, at low p the ground state 
only contains a dimers, whereas for higher p (> 0.3 for 
U = Q) both dimers are present. The resulting p{p) curve 
obtained from Had (not shown) exhibits an irregular pat¬ 
tern of steps of 2 or 1 particles. The two-component 
phase extends for values around U = 0, whereas for 
U Ji {U —Ji) only a {d) dimers are present [57]. 

The atom and doublon dimerizations, Na = 
T>i{blbi+i), and (bi+if) , approximately 

measure pa and pd (Fig. 3 (c)). Whereas the SF and PSF 
phases are LLs of almost hardcore single particles or pairs 
(for 0 = tt), in the PP phase both components coexist. 
Note the existence of many sharp jumps to states with 
approximately finite doublon density. 

The AHM is characterized by a density- and U - 
dependent drift of the quasimomentum distribution [30, 
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50], evident in Fig. 4 for the SF. Interestingly, whereas 
the SF presents an usual single-peaked distribution, the 
PP phase shows a multipeaked distribution. This is be¬ 
cause the a and d hop with a different Peierls phase, and 
hence the overall (a + d) momentum distribution presents 
multiple peaks. For PSF, single-particle correlations de¬ 
cay exponentially, and the peak structure blurs. 

Unit filling.-We consider now the special case of p = 1. 
Tuning 0 may induce a SF to MI transition [30, 50]. 
This is so even for 1/ <0 (Fig. 2 (c)) due to the effec¬ 
tive on-site repulsion discussed above. This Kosterlitz- 
Thouless transition may be determined as the curve at 
which K = 2 [58], which is consistent with the open¬ 
ing of the energy gap. As mentioned above, for 1/ <0 
an Ising-type transition occurs from SF to PSF. How¬ 
ever, in the vicinity of 0 ~ tt, the combination of anyonic 
statistics and 2BHCC results in a dimerized phase, D. 
The latter is characterized by a finite dimer-order pa- 
rameter Od = (Tl /2 - Tl/ 2 +i) with Ti = foj&i+i + H.c.. 
The Ising transition line is determined by the crossing of 
single- and two-particle excitations, and by a sharp peak 
in the fidelity susceptibility [59]. Interestingly, the phase 
diagram presents a SF-MI-D-PSF multicritical point. 

Probing the anyon-Hubbard model- A discussion on ex¬ 
perimental detection demands a grand-canonical analysis 
since local density arguments translate the dependence 
on the chemical potential, p, into a spatial dependence 
of the phases in the presence of an overall harmonic con¬ 
finement. Crucially, the PP phase has a finite extension 
in both (j) and p (Fig. 2 (d)). Hubbard models with corre¬ 
lated hopping (both for fermions [60, 61] and bosons [62]) 
may present a similar (but insulating) phase with coex¬ 
istent doublons and holons, which is, however, infinitely 
compressible, i.e. it is a single point as a function of ji. As 
a result these phases cannot be observed in trapped ex¬ 
periments, just leading to abrupt density jumps. On the 
contrary the novel PP phase, having a finite fi extension, 
may be experimentally observable by properly setting the 
central /i in the PP region (adjusting the particle num¬ 
ber). Moreover, the multipeaked momentum distribution 
of the PP phase can be easily probed in TOF measure¬ 
ments. Figure 4(c) shows that TOF images of A and the 
B may be employed to reveal the PP phase [32]. 

Summary. - We propose a Raman scheme that realizes 
the AHM using cold atoms in ID lattices, overcoming the 
main limitations of the proposal of Ref. [30]. This scheme 
results naturally in a 2BHCC, and allows for tailoring ef¬ 
fective on-site interactions without Feshbach resonances. 
The combination of anyonic character and 2BHCC leads 
to a rich ground state physics, including Mott insulators 
for attractive interactions, a dimerized gapped phase, 
pair superfiuids, and a novel two-component superfluid 
(partially paired phase) characterized by a large com¬ 
pressibility and a peculiar multipeaked momentum dis¬ 
tribution that may be easily probed experimentally. Fi¬ 
nally, we would like to mention that the Raman scheme 



FIG. 4: (a) Quasimomentum distribution n{k) for p = 1/2 
as function of I//t for 0 = O.Ott (L = 60 sites). Dashed lines 
denote phase transitions between PSF, PP and SF phases, 
(b) and (c) show the momentum distribution for U = 2t and 
U = 0 for the AHM (solid line) and the components A (dashed 
line) and B (dotted line) separately (for bosonic atoms). The 
latter are obtained by direct simulation of the two-component 
model with Raman-induced hopping (see Suppl. Material for 
a discussion of the fermionic case). 

may be extended to 2D lattices. Although for 2D lattices 
the boson-anyon mapping is lost, a 2D setup would open 
the exciting possibility of inducing a density-dependent 
magnetic field, contrary to the static field created in re¬ 
cent experiments [34, 63-65]. This possibility will be in¬ 
vestigated in a forthcoming work. 
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In this Supplementary Material we discuss in more extent the creation of density-dependent gauge fields using 
Raman-assisted hopping, and additional details about detection using the quasi-momentum distribution. 


I. CREATION OF DENSITY-DEPENDENT GAUGE FIELDS 

A. Two-component system 

We consider two hyperfine states, A and B, of a bosonic or 
fermionic atom. Both components are confined to the lowest 
band of a tilted ID optical lattice along x, of spacing D, and 
depth Vq = sEp, with Ep = l2mD the recoil energy. 
Without tilting there is a hopping rate J to nearest neighbors. 
The lattice tilting induces an energy shift A from site to site. 

We denote sls w{x — jD) the Wannier function maximally 
localized at site j. Due to the tilting, it is convenient to use 
Wannier-Stark states. For J <C A, the Wannier-Stark state 
centered at site I may be approximated as fji{x) w{x — 
ID) -1- ^{w{x — (/ + 1)^) “ ~ ~ 1)^))- The 3D on¬ 

site wavefunction at site j is ^j(r) = fjj{x)(p{y, z), where 
(p{y^z) is given by the strong transversal confinement. For 
simplicity we assume below ip{y^ z) w{y)w{z). 

On-site interactions between a and /S (a^/S =A,B) 

are characterized by the coupling constant Ua,p = 

I with aa,p the corresponding scatter¬ 

ing length. For a sufficiently deep lattice, the evaluation of 
the on-site interactions is simplified by means of the har¬ 
monic approximation: <F(r) , where 

I = Ds~^l^ j'K. Using this approximation we obtain: Ua,y — 

• Of course, for fermions only inter¬ 
species on-site interactions are possible. 


B. Laser transitions 

Before discussing Raman-assisted hopping, we comment 
on specific transitions in alkali atoms. We analyze first the 
case of bosons, which we illustrate with the example of ^^Rb 
(below we comment on other species). We choose A= \E = 
l^rriF = —1) and B= \E = 2, mi? = —2). Note that one 
could in principle swap the definition of A and B. The oppo¬ 
site choice would be however problematic, since the scattering 
lengths fulfill a{\E = l^mp = —1), \E = 2, mp = —2)) = 
a{\E = 2, mi? = —2), \E = 2, mi? = —2)) [1]. Hence swap¬ 
ping A and B would lead to aAA = clab, which would not 
allow to remove some spurious Raman transitions (see Sub¬ 
sec. IH). Note also that since we consider maximally stretched 
states, hyperfine-changing collisions are absent. 

We consider two lasers (see Fig. 1), one with linear po¬ 
larization (tt laser), and other with cr_ polarization (cr laser). 
We define as 6 the detuning of the tt laser from the |F = 
2, mi? = —2) ^ \E' = 2^m'p = —2) transition, which 
is (up to slight detunings of the order of kHz or smaller) the 


F' = 2 






FIG. 1: (Color online) Left: Transitions considered for the realization 
of the Raman-assisted hoppings, with tt polarized lasers (blue) and 
CT- polarized lasers (red); right: if x ^ 1 (see text), one may assume 
in a very good approximation that the tt lasers affect just the A state, 
whereas the a lasers affect just the B state by construction. 


same as that of the a laser from the \E = 1, m^ = —1) —> 

\E' = 2^m'p = —2) transition. Both lasers are red de¬ 

tuned. Note that due to selection rules \E = 2, mi? = —2) 
is dark with respect to the a laser. In contrast, the tt laser 

links as well \E = l,mi? = —1) to \E' = 2,m'^ = —1) and 

\E' = 1, m'p = — 1), with detunings, respectively, 6 + ojhfs 
and d^ujRFS ~^hfs^ where ujrfs and uj'^ps hyper¬ 

fine splittings of the ^5'i/2 and ^Pi /2 manifold, respectively. 
The fact that the tt laser affects \F = l,mi? = —1) may 
introduce spurious effects (see Subsec. ID). 

Fermions may be treated in a similar way. For the specific 
case of we choose A= \E = 9/2, mi? = —9/2) and 
B= \E = 7/2, mi? = —7/2). We may then employ a similar 
TT-cr laser configuration as in Fig. I. Since F = 9/2 lies below 
E = 7/2 the lasers must be now blue detuned. The above 
mentioned choice of A and B is better than the opposite, since 
it minimizes the spurious effects introduced the the coupling 
of \E = 7/2, mi? = —7/2) by the tt laser (see Subsec. ID). 


C. Raman-assisted hopping 


Direct hopping is prevented for a large-enough tilting, J <C 
A. Raman assisted hopping may be induced by the absorption 
of a photon of laser and the emission of a photon of laser 
Ln. The amplitude of the Raman-assisted hopping, Jnm^ from 
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Species 

(^HFS (GHz) 

^'hfs (GHz) 

r (MHz) 

X 

Na 

1.722 

0.189 

9.8 

0.21 


6.834 

0.814 

5.75 

0.059 

®®Rb 

3.036 

0.361 

5.75 

0.13 

Cs 

9.193 

1.167 

4.57 

0.046 


TABLE I: Values of uhfs, ^hfs ^ the linewidth L, and x (evaluated 
for = 0.5 GHz) for different atomic species. For Na and ®^Rb, A= 
\F — 2, mp — —2), and B= \F — l.rup — —1), and x is defined 
as in Eq. (4). For ®^Rb, A= \F — S.mp = —3), and B= \F = 

2,mF = -2), andx = 6 + 47 5+ppF^]- P”’’ 

Cs, A= \F = 4,mF = —4), and B= \F = 3,mF = —3), and 


X = 6 


[JL_5_ 

[48 S+ujhfS 


+ 


A_5_1 

16 S+ujhfS-^^'hfS \ ’ 


site j to site j + 1, may be evaluated following Ref. [2]: 

Jnra = ^ j ( 1 ) 


where is the difference between the laser 

wave vectors of the lasers, and Vnm is the amplitude of the 
Raman coupling (see Subsec. ID). For J ^ A\ 


/ 


dx'ipj+i{xye^^'^F^^.(^x) ~ 


. i 5 fe "’” D ( i + l / 2 ) 


J . 

2 ^ — sm 

A 


For s :§> 1 <C D): 

T ^ 

J n.m. — ^ 


-!>)•» W 




5k^^D 


2^i5k2^x 


2 A 


J sin 


5k^^D 


^i8k^^DU+ll2) 


( 2 ) 


. ( 3 ) 


Note that 5k'^'^ 7 ^ 0 is necessary to establish a significant 
hopping {2-A]. 


D. Four-laser scheme 

We employ the notation of the main text, to 

denote the states of neighboring sites, where may be 

0 (empty site), A or B (singly occupied site with an A or B 
component), and AA, AB and BB for doubly-occupied sites. 


1. Bosonic case 

Let us first study the bosonic case. We consider the Raman 
scheme of Fig. 1 of the main text, where ^ 2,3 are a lasers, 
and 1 / 1^4 are tt lasers. The lasers are characterized by their 
Rabi frequency proportional to the square root of 

their intensity, their frequency ujj , and their wave vector kj. 
As mentioned above, we illustrate the case of bosonic species 


with the example of ^^Rb, with A= \F = 2^mF = —2), 
and B= \F = l^rriF = —1). Let us consider in detail the 
different Raman-assisted hops we would like to induce (the 
Clebsch-Gordan factors may variate for other alkali species): 

• (i) (A, 0) ^ (0, A). This transition leads to an energy 
shift AE = —A. For co ’2 — ^3 = — ^4 = —A, the 

transition is realized by either Li and I/ 4 , or by L 2 and 
I/ 3 , with amplitude V 23 = \ with 


X = 6 


"1 ^ 

_4 ^ + lchfs 


^_ 5_ _ 

12 ^ + uhfs — ^hfs 


(4) 


• (ii) {A, A) ^ (O^AB). This transition has AE = 
Uab —As., and involves absorbing a a photon, and emit¬ 
ting a TT one. We assume uj 2 —^4 = Uab —A-\-U, where 
(7 is a small detuning (within the Raman resonance 
width). As discussed below this tunneling plays the role 
of effective on-site interactions in the final model. The 
transition is quasi-resonantly provided by L 2 and I/ 4 , 
with amplitude ^24 = 

• (m){AB,0) {A, A). For this case, AE = Uab —As, 

and is performed by absorbing a tt photon, and then 
emitting a a one. For cci — CC 3 = —Uab — A — U, the 
transition is quasi-resonantly provided by Li and I/ 3 , 
with Raman amplitude V 13 = ^ . 

• (iv) {AB,A) ^ {A,AB). For this transition AE = 
—A, and can be only realized using tt lasers. For 
uji — UJ 4 = —A, the transition is realized by ^ 2 , 3 , with 

amplitude 1^4 = ^ . 

|f^l||f^4| _ |f^2||f^3| _ l^lll^sl _ |f^2||f^4| _ 02 

f 2 i = |r 2 i|e“®'^, and $7^=2,3,4 = one obtains V23 = 

^(1 + V24 = Vi3 = and 

Li 4 = The term proportional to x in (i) is not 

discussed in the main text, for simplicity of the exposition. 
It results in an additional term in the final Hamiltonian, as 
shown in Subsec. IE. However, x can be kept small. Specific 
examples for different species for ^ = 0.2 GHz are shown in 
Table I [5-8]. We show below that for small x the spurious 
term does not modify the conclusions of the main text. 


2. Fermionic case 


For fermions it is better to choose A as the state with larger 
F, e.g. A= \F = 9/2,-9/2) and B= \F = 7/2,-7/2) 
for In the scheme of Fig. 1 of the main text, we choose 
1 / 2,3 as TT lasers, and I/i ^4 as a lasers. We may proceed as for 
bosons, but choosing fti = — Due to the inverted 

choice of A and B, we obtain: V 23 = X’ ^24 = 

Vi 3 = V 2 ^e-^^, and Vu = 2 ^( 6 “^^ -h x)- For 

1 

with ujhfs = 


32 


49 


^ [243 SF^hfs 243 SF^hfs—^'hfs 

1.286 GHz, OiJFFS — 0-155 GHz (the linewidth is F = 5.96 





3 


MHz) [9]. For S = 0.2 GHz, x = 0.072. Note that inverting 
the choice of A and B, would lead to Vnm of the same form as 
in the bosonic case but with x four times larger (x = 0.28 for 
S = 0.2GHz). The latter justifies our choice of A and B. 


E. Effective Hamiltonian 

We consider that only processes (i) to (iv) are resonant, and 
postpone to the Subsec. IH the discussion about spurious pro¬ 
cesses. We denote as cj the bosonic operator corresponding 
to the Fock-state manifold {10), |1) = \A), | 2 ) = \AB)}. Let 
us consider ki 2 = key, and ks 4 = kex, and k = ir/D. 


1. Bosonic case 

For the bosonic case, the Raman-assisted hopping of the 
previous subsection result in the effective Hamiltonian: 



FIG. 2: Equation of state p = p{p) of the AHM {U = 0, 0 = tt, L = 
48) for different values of x for the spurious correlated tunneling in 
(a) Eq. (6) and (b) Eq. (7). Small arrows depict the kinks in the p{p) 
curve for X = 0.2 characterizing the PP-SE transition. 


H = 




+ - nj+i)cj +R.C. 


U 


- 1 ). 


( 5 ) 


with Uj = YjCj, and t = J- Typical values of the 

Raman-assisted hopping rate, t, are of the order of few tens 
of Hz. Note that the factor (—1)-^ , which results from the 
X projection of Sk, may be easily eliminated by redefining 
the bosonic operators in the form: b^i = C 41 , 64/+1 = C 4 ^+i, 
hi +2 = -C 4 Z+ 2 , ^ 4 Z +3 = -C 4 /+ 3 , with I an integer. In this 
way we obtain: 


j 

+ xE - «i)(l - + H.c.' 

j 

j 


F. Effect of a finite x > 0 

For X = 0, Hamiltonians ( 6 ) and (7) reduce to Eq. (1) in 
the main text, and hence to the anyon-Hubbard model (AHM). 
X > 0 results in spurious collisionally-assisted hopping terms. 
However, as commented above, x can be kept small. 

Figure 2 (a) shows the effect of the spurious correlated tun¬ 
neling in Eq. ( 6 ) which basically effects the low-density part 
p < 1 and results in a small shift of the SF-PP transition (in¬ 
dicated by small arrows) to lower p. Interestingly a PP phase 
may be observed also for p > 1 for sufficiently large x- 

Figure 2 (b) shows the effect of the spurious correlated tun¬ 
neling in Eq. (7). The effect is negligible for small lattice fill¬ 
ings p < 1. Due to its p^ -dependence even for much larger 
values of x as discussed in the previous subsection the curves 
of Fig. 2 (b) basically overlap with the pure AHM-case. For 
densities p > 1 the phase boundaries are slightly modified. In 
particular the size of the Mott-insulator region increases and 
the PP phase may now extend to a small region above unit 
filling (see arrows in Fig. 2 (b)). 

Note that, interestingly, the choice of a larger detuning 5 
and thus larger values of x actually favors the observation of 
the PP phase in experiments. 


2. Fermionic case 

Proceeding in exactly the same way for the fermionic case, 
we obtain 

3 

u . 


G. Comparison to the model of Keilmann et al. 

Note that (i) and (iv) are energetically degenerated, but 
(apart from the small spurious effect discussed above) they 
may be addressed with different lasers due to selection rules. 
This point constitutes the major drawback of the proposal of 
Ref. [10]. In that proposal, a single component (A) was con¬ 
sidered, and process (iv) was of the form (AA,A)^(A,AA), 
which cannot be resolved from process (i). As a result, both 
(^) 1/2 and I/ 3 , and Li and I /4 address both (i) and (iv), preventing 

the realization of the desired density-dependent Peierls phase. 
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Species 

claaIclb 

^abI^b 

\5U/h\ (kHz) 

Na 

54.53 

64.5 

0.72 


100.36 

99.08 

0.23 

®®Rb 

-469 

-386 

1.40 

Cs 

-1980 

2436 

8.89 


TABLE II: Values of a a a, (iab, and 5U /h for different bosonic 
species, for s = 20 and D = 532nm (except for where we 
employ s = 25 and D = 266nm) [14]. For Na and ®^Rb, A= 
\F = 2,mF — —2) and B= \F = l^mp — —1). For ®^Rb, 
A= \F = 3,mF = —3), and B= \F = 2,mF = —2). For Cs, 
A= \F = — —4),andB= \F = 3,mF = —3). Note that 

negative scattering lengths are of no concern in our set-up. 


The two processes may be just discerned by considering a very 
small detuning 6 A. Since Uaa and A are of the or¬ 

der of kHz, 6 < 100 Hz is necessary. However F is typically 
of the order of few MHz (see e.g. Table I), rendering unfeasi¬ 
ble the proposal of Ref. [10]. 

H. Spurious processes 

Spurious Raman processes must be kept out of resonance: 

• (v) (A,0) ^ (0,B): AE = -A 

• (vi) (A,A) ^ ( 0 ,AA): AE = -Af Uaa 

• (vii) (AA,0) ^ (A,A): AE = -A- Uaa 

• (viii) (AB,A) ^ (B,AA): AE = -Af SU, with SU = 
{Uaa - Uab) 

• (ix) (AA,B) -> (A,AB): AE = -A- SU 

Note that (vi), (vii), (viii) and (ix) are not possible for 
fermionic atoms. Process (v) is just possible with J 24 or 
J 13 . But these laser combinations are (quasi-)resonant with 
—A ± Uab- For a sufficiently large Uab process (v) is hence 
far from resonance with either J 24 or J 13 . As an illustra¬ 
tion we consider the case of atoms, with aAB — ITOa^, 
with as the Bohr radius, and a lattice with D = 532 nm, 
and s = 20, Uab/^ ^ 7 kHz. Note that typical values for 
the width of the Raman resonance, W < 50 Hz [2]. Hence 
Uab ^ W, and hence process (v) may be safely ignored. 

For bosons, (vi), (vii), (viii) and (ix) must be also consid¬ 
ered. These processes can be only realized using J 23 . To ne¬ 
glect the (vi) and (viii) processes one needs Uaa ^ W. This 
is the case in experiments, where typically Uaa ^ 1 to few 
10s kHz. In contrast, to avoid (viii) and (ix) one must demand 
SU ^ W. The latter condition is more strict, but may be at¬ 
tained as well in experiments. Examples for different bosonic 
species are illustrated in Table II [14, 15]. 


I. Initial preparation 

If the spurious processes can be neglected, only on-site 
states 0, A, and AB may be reached by Raman-assisted hop¬ 
ping. However, B, AA, BB, and even sites with occupation 
higher than two, may initially exist. Sites with an occupation 
larger than two are eventually destroyed by three-body pro¬ 
cesses, and are hence of no major concern. This is however 
not the case of B, AA and BB, which should be avoided in 
the initial preparation, since they would induce immobile 
spurious defects. A straightforward way to avoid immobile 
defects is to consider a low-density gas of maximally singly 
occupied sites by A particles (doubly occupied AA sites may 
be eliminated using photo association pulses [13]). One may 
then compress this dilute gas increasing the density until 
reaching the desired density. Since A can just move into 
AB, doubly occupied sites will be guaranteed to have an AB 
occupation, and hence the on-site states remain always within 
the {0, A, AB} manifold. 


II. QUASI-MOMENTUM DISTRIBUTION 

As discussed above, the AHM may be realized experimen¬ 
tally with both fermionic and bosonic species. While many 
observables, most relevantly the density, do not depend on the 
particular choice of constituents, others, in particular the mo¬ 
mentum distribution, are not gauge invariant. 

As shown in Fig. 4(b) and (c) of the main text, bosonic 
species closely mimic the quasi-momentum distribution of the 
effective model. As mentioned above, the oscillating factors 
( —1)-^ resulting from the necessary x projection of ^k, may 
be easily eliminated by a gauge transformation. However, this 
transformation must be undone in order to recover the mo¬ 
mentum distribution of components A and B. As a result, A 
and B components exhibit a doubled unit-cell, and the quasi¬ 
momentum distribution n{k) is consequently trivially shifted 
and duplicated. For clarity we only depict half the Brillouin 
zone for A and B in Fig. 4 of the main text. 

It is interesting to depict as well the momentum distribu¬ 
tion for the fermionic case. Figures 3 (a) and (b), which de¬ 
pict n{k) neglecting the ( —1)-^ factors, show a single Fermi- 
surface in the SF, two for the PP, and a blurred n{k) for the 
PSF. The oscillating (—1)-^ factors obscure this simple inter¬ 
pretation (as shown in Figs, (c) and (d)). However PSF, PP 
and SF phases may be still clearly discriminated from each 
other by their characteristic momentum distribution function. 


[1] The s -wave collision between \F = l^mp = —1) and \F = 
2,mF = —2) for homonuclear collisions can occur only on 


the triplet curve. Hence, the scattering length at zero field is 
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FIG. 3: Quasi-momentum distribution n{k) of the A component for the case of fermionic atoms as function ofU jt {(j) — ti, p — 1/2, L = 60 
sites), (c) and (d) show n{k) including the oscillating factors (—1)^ created by the Raman-assisted hopping; (a) and (b) show the same data, 
but excluding the oscillating terms, (b) and (d) show n(k) for t/ = 2 (solid line), f/ = 0 (dashed line) and U — —2 (dotted line). In (a) and 
(c) the dashed lines denote the phase transitions between PSF, PP and SF phases. 
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